Resonant Enhancement of Inelastic Light Scattering in Strongly Correlated Materials 
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We use dynamical mean field theory to find an exact solution for inelastic light scattering in 
strongly correlated materials such as those near a quantum-critical metal-insulator transition. We 
evaluate the results for q = (Raman) scattering and find that resonant effects can be quite 
large, and yield a triple resonance, a significant enhancement of nonresonant scattering peaks, a 
joint resonance of both peaks when the incident photon frequency is on the order of U, and the 
appearance of an isosbestic point in all symmetry channels for an intermediate range of incident 
photon frequencies. 
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Inelastic light scattering is a powerful tool to unravel 
the nature of elementary excitations in a wide variety 
of materials [1], ranging from Kondo insulators [2, 3], 
to high temperature superconductors [4, 5], to colossal 
magnetoresistance materials [6]. The experimental ef- 
forts have grown tremendously with the availability of 
third generation light source facilities and improvements 
in CCD detectors. These efforts have been brought to 
bear on strongly correlated materials to examine the el- 
ementary excitations of insulators and metals and how 
they evolve as the correlations are made to change via 
doping, for example. 

One of the most studied areas is how the scattering 
cross section resonates with the incoming light frequency. 
It is widely believed that by tuning the incident pho- 
ton frequency, features of the non-resonant spectra can 
be magnified by orders of magnitude; that is, the reso- 
nance serves as a bootstrap to raise the intensity of the 
non-resonant signal. However, a full, consistent theory is 
lacking [7, 8] . Non-resonant scattering is derivable from a 
two-particle correlation function which can be treated by 
a variety of techniques, yet the resonant and mixed con- 
tributions involve higher particle correlations and are dif- 
ficult to treat theoretically due to multiple-particle vertex 
renormalizations. Most of the approaches to light scat- 
tering in insulators examine the Loudon-Fleury model [9] 
which is appropriate for off-resonant conditions for the 
scattering of light off of spin excitations, for example. In 
the strong-coupling regime, a perturbative approach has 
been used to illustrate a number of important features 
of electronic resonant scattering processes [8, 10]. The 
nonresonant case has also been examined, and an exact 
solution for correlated systems (in large spatial dimen- 
sions) is available for both the Falicov-Kimball [11] and 



Hubbard [12] models. Here we concentrate on an exact 
solution of the full problem for the Falicov-Kimball model 
including all non-resonant, resonant and mixed scatter- 
ing channels. 

For an electronic system with nearest-neighbor hop- 
ping, the interaction with a weak external transverse elec- 
tromagnetic field A is described by [8] 

2 

H int = ~j-A+^J2 A ^«pA , (1) 
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where 

ja = Mk)ct{k + q/2)c a (k - q/2), (2) 

k 

is the current operator, v a (k) = de(k) / dk a is the Fermi 
velocity, and 

™ = E |r^r c ^ fc + 9/ 2 ) c -( fc - 9/ 2 ) ( 3 ) 

k a ^ 

is the stress tensor operator. The inelastic light- 
scattering cross section becomes (fl — u>i — ujf, q = 
ki — kf is the transferred photon frequency and momen- 
tum, respectively): 

R(fl) = R N (£l) + Rm(Q) + R R (n), (4) 

where the nonresonant contribution is 

R N (fl) = 27rg 2 (k z )g 2 (k f ) 

x V ~ 7i/ 7/i S(e f -Si- 0), (5) 
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the mixed contribution is 
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and the resonant contribution is 
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Here 0Ji(f) and fcj(^) denote the energy and momentum of 
the initial (final) states of the photons, £ i( ^) are the eigen- 
values corresponding to the eigenstates that describe the 
"electronic matter", and g(k) = (he 2 /VuJk) 1 ^ 2 is the 
"scattering strength" with uou = c\k\. We have intro- 
duced the following symbols 



7 = H e aTa/3e^, 
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with the notation Oij — (i\0\f) for the matrix elements 
of an operator O, Z the partition function, and e 1 '? are 
the incident and scattered light polarization vectors, re- 
spectively. We concentrate on the light scattering re- 
sponse function x(^)i which is related to the cross sec- 
tion, but with a Bose statistical factor removed: 



(9) 



1 - exp(-/3fi) 

Inelastic light scattering examines charge excitations of 
different symmetries by employing polarizers on both the 
incident and scattered light. The A ig symmetry has the 
full symmetry of the lattice and is primarily measured 
by taking the initial and final polarizations to be e* = 
= (1,1,1,...). The Bi g symmetry involves crossed 
polarizers: e 4 = (1,1,1,...) and = (—1,1,-1,1,...); 
while the £?2 g symmetry is rotated by 45 degrees, with 
e J = (1,0,1,0,...) and = (0,1,0,1,...). While a sym- 
metry analysis can be employed for all momentum trans- 
fers q, for Raman (q = 0) scattering, it is easy to show 
that for a system with only nearest neighbor hopping and 
in the limit of large dimensions, the A\ g sector has contri- 
butions from nonresonant, mixed, and resonant scatter- 
ing, the Bi g sector has contributions from nonresonant 



and resonant scattering only, and the i?2 g sector is purely 
resonant [11]. The symmetry analysis would be substan- 
tially different for lower dimensions but is not currently 
tractable. A full analysis for all q will be presented else- 
where and therefore for the remainder of the paper we 
focus on Raman scattering (q = 0) only. 

Normally the matrix elements defined in Eq. (8) can- 
not be easily determined for a many-body system in the 
thermodynamic limit. Instead, the light scattering cross 
section expressions must be evaluated by first consider- 
ing the relevant multi-time correlation functions on the 
imaginary time axis, then Fourier transforming to a Mat- 
subara frequency representation, and finally making an 
analytic continuation from the imaginary to the real fre- 
quency axis. In the case of nonresonant scattering, the 
expressions to be analytically continued depend on only 
one frequency; for mixed scattering they depend on two 
frequencies, and for resonant scattering, they depend on 
three. The analytic continuation procedure for the mixed 
and resonant Raman scattering is complicated, because 
it requires a multistep procedure, where first the trans- 
ferred frequency is continued to the real axis, then the 
individual initial and final frequencies are continued to 
the real axis. In addition to the analytic continuation, 
we also must evaluate the dressed multi-time correla- 
tion functions. There are renormalizations associated 
with two-particle "ladder-like" summations for a num- 
ber of the relevant diagrams, but the symmetry of the 
velocity operator, and of the relevant multi-particle ver- 
tex functions (which are local in the large-dimensional 
limit) imply that there are no parquet-like summations, 
nor are there any three- or four-particle vertex renormal- 
izations [13]. Since the two-particle vertex function for 
the Falicov-Kimball model is already known [14], the full 
Raman scattering problem can be solved via a straight- 
forward but tedious procedure. The final formulas are 
cumbersome and will be presented elsewhere. 

The Falicov-Kimball model Hamiltonian satisfies [15] 

H = -TTnY.^ + c h) + uJ2 c t c * w > ( 10 ) 

(ij) » 

where c\ (ci) create (annihilate) a conduction electron at 
site i, u>i is a classical variable (representing the localized 
electron number at site i) that equals or 1, t* is a renor- 
malized hopping matrix that is nonzero between nearest 
neighbors on a hypercubic lattice in d-dimensions (and we 
take the limit d — > oo [16]), and U is the local screened 
Coulomb interaction between conduction and localized 
electrons. This model can be solved exactly by using dy- 
namical mean field theory, as described by Brandt and 
Mielsch [17] and summarized in review articles [18]. 

We concentrate on the case with U = 2 here, which 
is just on the insulating side of the metal-insulator tran- 
sition at half filling (p e — (wi) = 1/2). This was the 
regime where the nonresonant response showed a num- 
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ber of interesting properties for both Raman [11] and 
inelastic x-ray scattering [19]. The Stokes Raman re- 
sponse function is plotted in Fig. 1 at T = 0.5 for 9 
different incident photon frequencies uii ranging from 0.5 
to 4.5 in steps of 0.5. Since the transferred energy can be 
no larger than the incident photon energy all scattering 
curves run from zero up to u>i. The first thing to note 
in Fig. 1 is the large nearly vertical line that occurs as 
— ► Wi. This is the triple resonance [10], which yields a 
strong enhancement to the Raman scattering when the 
energy of the scattered photon approaches zero. In the 
Loudon-Fleury regime [9] , where the incident photon en- 
ergy is much larger than the electronic excitation ener- 
gies, we see that the full response is essentially that of the 
nonresonant response [11] plus the triple- resonance peak. 
As the incident photon energy is reduced, the behavior 
becomes much more complex. Gcncrically we can iden- 
tify a number of low-energy and higher-energy resonant 
enhancements to the scattering. 
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FIG. 1: Raman response function for different channels. We 
take U — 2, T = 0.5, and choose w< = 0.5 — 4.5 in steps of 0.5 
(the different line thicknesses correspond to different w»'s). 

One interesting feature of the response function, seen 
in experiments on correlated materials [2, 3], and seen in 
theoretical calculations of the nonresonant response [11, 
12], is that at low energy there is an isosbestic point, 
where the B\ s response function is essentially indepen- 
dent of temperature at a particular frequency 17 « U/2. 
Below that frequency the response decreases as T is low- 
ered, and above it increases. The isosbestic behavior 
must survive in the Loudon-Fleury regime, because the 
isosbestic point is at low energy, and the low-energy re- 
sponse is negligible in the resonant and mixed contribu- 
tions. But what happens when Ui w Ul Here we expect 
interesting effects to occur, because the incident photon 



energy is the right size to cause transitions from the lower 
to upper Hubbard bands of the correlated insulator. In- 
deed, we find interesting results in this regime (Fig. 2). 
At low temperature [T < 0.7), a symmetry-dependent 
isosbestic point appears at a transferred frequency of 0.7- 
0.9 and is seen in all channels at low enough T, even the 
Ai g and Bi g channels, which have no isosbestic point 
in the nonresonant regime. Hence the inclusion of res- 
onant and mixed terms provides theoretical support for 
the generic presence of a low-temperature isosbestic point 
in correlated systems. 
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FIG. 2: Raman response function for U = 2 and fi = 2 for 
different channels at T = 1, T = 0.5, T = 0.2, and T = 0.05. 
The temperature decreases as the lines are made thicker. 

Finally, we present results of what the resonant pro- 
file of the scattering looks like by fixing the transferred 
frequency and varying the incident photon energy. We 
expect that there will be a resonant peak in the response, 
and indeed this is so, although in some cases the triple 
resonance overwhelms the presence of the peak. Note 
that in high-temperature superconductors, in addition 
to the expected resonance that occurs when the incident 
photon frequency is close to the transferred frequency, 
another resonance occurs, where the low-energy peak is 
strongly enhanced when u>i ss U [4, 5]. We show this 
regime in Figs. 3 (O = 2.0) and 4 (Q. = 0.5). In Fig. 3 
we see a moderately broad peak centered at Wj 10-20% 
higher than U. The enhancement of the charge-transfer 
peak in this regime can easily be an order of magnitude 
over the nonresonant response. In Fig. 4, we see a simi- 
lar resonant feature when the incident photon frequency 
is slightly larger than Q = 0.5 (arising from the triple- 
resonance effect), but a second less prominent series of 
broad peaks occurs when u>i w U indicating that the 
low-energy and charge transfer peaks are resonating to- 
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gethcr when Ui w U . Hence the behavior observed in 
the high-temperature superconductors [4, 5] is likely to 
be seen in many other correlated insulators. 
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FIG. 3: Raman response function for U = 2 and Q = 2 for 
different channels at T = 1, T = 0.5, T = 0.2, and T = 0.05. 
The temperature decreases as the lines are made thicker. 
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FIG. 4: Raman response function for U = 2 and Q = 0.5 for 
different channels at T = 1, T = 0.5, T = 0.2, and T = 0.05. 
The temperature decreases as the lines are made thicker. 

In conclusion, we have shown a number of interesting 
resonant features in theoretical calculations of electronic 
Raman scattering. These features include the triple res- 
onance, the resonant enhancement of nonresonant peaks, 
the appearance of isosbestic points, and the joint res- 
onance of low-energy and charge-transfer peaks when 



LUi w U . It will be interesting to see whether these pre- 
dictions can be seen in future experiments on correlated 
systems. 
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